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Gauge Coupling Beta Functions in the Standard Model to Three Loops 
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In this Letter we compute the three-loop corrections to the beta functions of the three gauge 
couplings in the Standard Model of particle physics using the minimal subtraction scheme and 
taking into account Yukawa and Higgs self-couplings. 
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Renormalization group functions are fundamental 
quantities of each quantum field theory and play an im- 
portant role in various aspects. Besides controlling the 
energy dependence of parameters and fields they are also 
crucial for the resummation of large logarithms. Further- 
more, renormalization group functions are important for 
the development of grand unified theories and the extrap- 
olation of low-energy precision data to high energies, not 
accessible by collider experiments. 

As far as the strong interaction part of the Standard 
Model is concerned the corresponding gauge coup ling 
beta function is known up to four- loop order [lMlOj. At 
three-loop level also the corrections involving two strong 
and one top quark Yukawa coupling have been com- 
puted [HI. On the other hand, for the SU(2) L x U(l) Y 
part only one- [l], 0, H3| and two-loop |13l - [l6j corrections 
are available since about 30 years. (Two-loop correc- 
tions are also known for the Yukawa |16l - tl8| and Higgs 
boson self-coupling 0, Ql| [20| beta function, see also 
Ref. [H]].) For a general theory based on a simple gauge 
group the three-loop corrections to the gauge coupling 
beta function have been calculated in Ref. |22|. In this 
Letter we provide results for the three-loop gauge cou- 
pling beta functions taking into account all sectors of the 
Standard Model, i.e., the gauge, Yukawa and Higgs boson 
self-couplings. 

Let us in a first step define the beta functions. We 
denote the three gauge couplings by a,\, a 2 and a^ and 
adopt a 5(7(5)-like normalization with 
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where a is the fine structure constant, 9w the weak 
mixing angle and a s the strong coupling. In our cal- 
culation we consider in addition to the gauge couplings 
also the third-generation Yukawa couplings [34j a^ = a t , 
as = off, and a^ = a T , and the Higgs boson self-coupling 
a 7 = \/(Att). a x = am x / (2 sin 2 9 W M W ) (x = t,b,r) 
where m x and M\y are the fermion and W boson mass, 
respectively, and — A($t$) 2 j s the part of the Lagrange 
density describing the quartic Higgs self interaction. 

The functions fti are obtained from the renormalization 
constants of the corresponding gauge couplings that are 
defined as g^ aYO = [fZ gi gi where on = g 2 /(47r). Exploit- 



ing the fact that the <?k arc are /i-independent and taking 
into account that Z gi may depend on all seven couplings 
leads to the following formula 
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where e = (4 — d)/2 is the regulator of Dimensional Reg- 
ularization with d being the space-time dimension used 
for the evaluation of the momentum integrals and the 
dependence of on the renormalization scale (i is sup- 
pressed. From Eq. ([2]) it is clear that the renormalization 
constants Z gi (i = 1,2,3) have to be computed up to 
three-loop order. 

In the modified minimal subtraction (MS) renormal- 
ization scheme the perturbative expansion of the gauge 
coupling beta functions can be written as 
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In this Letter we evaluate the three-loop terms (coeffi- 
cients Cijk) only for the gauge couplings (i.e. i = 1,2,3). 
For our calculation the beta functions for the Yukawa 
couplings are needed to the one-loop order and the tree- 
level expression [first term in Eq. ((5}] is sufficient for f3\. 

In the MS scheme the beta functions are mass inde- 
pendent that allows us to use the Standard Model in the 
unbroken phase as a framework for our calculation. By 
construction less vertices have to be considered than after 
spontaneous symmetry breaking. The electroweak gauge 
bosons are denoted by W and B corresponding to the 
SU(2)l x U{l)y gauge groups, respectively. 

In principle each vertex containing the gauge coupling 



at tree level can be used in order to obtain Z gi via 
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where Z vert stands for the renormalization constant of 
the vertex and Zk,yrf for the wave function renormaliza- 
tion constant; k runs over all external particles. We have 
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FIG. 1: Sample Feynman diagrams contributing to the 
W1W2W3 vertex [(a), (b) and (c)] and the ghost-W^ vertex 
[(d) and (e)]. Dashed, dotted, solid and wavy lines corre- 
spond to Higgs boson, ghost, fermion and W boson lines, 
respectively. Sample diagrams for gluonic Green's functions 
are obtained for (a), (b), (d) and (e) by replacing the W lines 
by gluon lines. 



computed Z gi using both the ghost-gluon and the three- 
gluon vertex and Z g2 has been evaluated with the help 
of the ghost-W-3, the W1W2W3 and the 4> + 4>~W 3 vertex 
where (jy^ is the charged component of the Higgs dou- 
blet corresponding to the Goldstone boson in the broken 
phase and W%, W2 and W3 are the components of the 
W boson. Both for Z g2 and Z g3 the different ways lead 
to the same result that constitutes a strong check for 
the correctness of the final result. Z gi is obtained from 
vertices containing the B boson. Because of the Ward 
identity there is a cancellation of Z V ert and the factors 
\/ Zk,w{ other than the one corresponding to the B bo- 
son. Thus Z gi is solely computed from the wave func- 
tion renormalization constant of the B boson. Several 
three-loop sample diagrams contributing to the consid- 
ered three-point functions are shown in Fig. [TJ Because 
of the fact that the /3j do not depend on any kinemat- 
ical scale we evaluate the vertex functions in the limit 
where one external momentum is set to zero. In this way 
all loop-integrals are mapped to massless two-point func- 
tions that up to three loops can be computed with the 
help of MINCER [HI. 

An important issue in the present calculation is the 
treatment of 75 within dimensional regularization. Non- 
trivial contributions may arise if in the course of the cal- 
culation two fermion traces occur where both of them 
contain an odd number of 75 matrices and four or more 



7 matrices. It is straightforward to see that the Green's 
functions that we have chosen for calculating the beta 
functions contain at most one-loop triangle sub-diagrams 
(see, e.g., Fig. HJa)). This could potentially lead to con- 
tributions where a careful treatment of 75 is required. In 
our case, however, all these contributions vanish identi- 
cally due to anomaly cancellations [see, e.g., Ref. [zij ] 
since we always sum over all fermions of one generation. 
This has also been checked by an explicit calculation (2|| . 

Our calculation is based on a high level of automation 
in order to avoid errors due to manual interaction. As 
a core of our set-up we use the well-tested chain of pro- 
grams that work hand-in-hand: QGRAF [26[ generates all 
contributing Feynman diagrams. The output is passed 
via q2e to exp [27|, l28( that generates the FORM [29j code. 
The latter is processed by MINCER [23J that computes the 
Feynman integrals and outputs the e expansion of the 
result. A serious problem that had to be overcome in the 
course of the present calculation is the enormous number 
of diagrams that contribute to the individual renormal- 
ization factors. For example, in the case of the W\ W2W3 
vertex one has about 380 000 diagrams for the vertex it- 
self and 60 000 for the W3 self energy with similar num- 
bers for the other Green's functions. Thus, in total more 
than one million diagrams have to be computed. In order 
to handle such an amount of diagrams we have decided 
to automatically split the original problems into blocks 
containing only of the order of 1000 Feynman diagrams. 
In total we end up with a CPU time of about 100 days on 
a single core. Since our calculation is highly paralleliz- 
able, the final wall-clock time is about one day on 100 
cores. 

Note that at three-loop level all sectors of the Stan- 
dard Model contribute to the (3 functions of the gauge 
couplings. Thus the huge number of diagrams mainly re- 
sults from the numerous interaction vertices in the Stan- 
dard Model. In the set-up described above the strong 
interaction part of the Standard Model has already been 
used for a variety of calculations. We have added the 
electroweak part by establishing an interface to a model 
file of FeynArts [30j. It generates all Standard Model 
Feynman rules in a format readable by QGRAF and q2e. 
The model file for FeynArts has been generated using 
FeynRules [3l|. 

In the following we present results for j3\ , P2 and up 
to three loops. In order to keep the expressions compact 
we set the Yukawa couplings ab and a T to zero. Also, 
we use unit Cabibbo-Kobayashi-Maskawa matrix in our 
calculations. The complete result with ab and a T kept 
non-zero can be found elsewhere [25| . Omitting the /1 
dependence in the arguments of ctj (and taking the limit 
e — >• 0) we obtain 
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where rig labels the number of generations and n t the 
number of heavy up-type quarks. In the Standard Model 
we have no = 3 and n t — 1. The result including a fourth 
generation 35] (with Yukawa coupling ar) is obtained 
with the replacements (a t n t ) n — > (a t + at')™ (n = 1,2) 
and a 2 nt — > a 2 + a 2 ,. It is interesting to note that, 
although the two-loop expression is A-independent, the 
three-loop term of /3\ and 02 contains both linear and 
quadratic terms in A. The latter arise from diagrams as 
those in Fig. QJc). 

There are several checks on the correctness of our re- 
sult. The interface between FeynArts and QGRAF/q2e 
has been checked by evaluating several one- and two-loop 
results that are known in the literature. We have even 
considered quantities within the Minimal Supersymmet- 
ric Standard Model, like the relation between the squark 
masses within one generation, which is quite involved in 
ease clcctroweak interactions are kept non-zero. We have 
furthermore reproduced the two-loop results in the litera- 
ture both for th e gauge [l3T[l6l | and the Yukawa coupling 
beta functions [16l4l8| that constitutes a strong check 
on the correctness of the Feynman rules and the general 
set-up. This is also true for the three-loop corrections 
to the BBB vertex which we have verified to be zero. 
Also the three-loop corrections of order a^ctt to 03 from 
Ref. [ll| have been reproduced. All our calculations have 
been performed for three general gauge parameters, one 
for each gauge group. Whereas the renormalization con- 
stants for the vertices and wave functions still depend on 



the gauge parameters the Z gi are independent leading to 
gauge parameter independent results for . In a further 
check we have ensured that in the vertex diagrams no 
infra-red divergence is introduced although one external 
momentum is set to zero |36J ■ This is done by assigning a 
non- vanishing mass m to internal lines (a common mass 
for all particles is sufficient) and performing an asymp- 
totic expansion 32] in the limit q 2 3> m 2 where q is the 
non- vanishing external momentum of the vertex diagram. 
Asymptotic expansion for this limit is automated in the 
program exp and thus can be performed with the set-up 
described above. It is sufficient to restrict to the leading 
term in m 2 jq 2 and check that no ln(m 2 /^ 2 ) terms appear 
in the final result. Nevertheless the calculation becomes 
significantly more complex; some diagrams develop up to 
35 sub-diagrams when applying the rules of asymptotic 
expansion. With this method we have explicitly checked 
that the WiW 2 W 3 and three-gluon vertex are free from 
infra-red divergences. Since the results agree with the 
ones obtained from the other vertices also the latter are 
infra-red safe. 

In order to estimate the numerical effect of the new 
terms we use the experimentally measured values for 
ai(Mz) (i = 1,2,3) and run to a higher scale \i using 
two or three loops for the gauge and two loops for top 
Yukawa beta functions. For /.i = 2 TeV the relative dif- 
ference between two and three loops amounts to 0.003%, 
0.010% and 0.005% for cti, a 2 and a 3 , respectively. In 
the case of a 3 the shift is significantly smaller than the 
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one introduced by the experimental uncertainty at the 
Z boson mass scale. On the other hand, for a.\ and 
the three-loop effect is of the same order of magnitude 
compared to the experimental uncertainty. Similar con- 
clusions hold for /i = 10 16 GeV where a relative shift be- 
tween two- and three-loop running amounts to 0.012%, 
0.027% and 0.004% for ai, ai and «3, respectively. 

To conclude, in this Letter we have computed the 
three-loop corrections to the gauge coupling beta func- 
tions in the Standard Model. Whereas Yukawa correc- 
tions are already present at two-loop order the Higgs bo- 
son self-coupling appears for the first time at three loops. 
The numerical effect of the new terms is small, however, 
in the case of u\ and oli it is comparable to the experi- 
mental uncertainties. 

The method used for the current calculation can also 



be applied to the calculation of the beta functions of the 
Yukawa and Higgs boson self-coupling. However, in the 
case of the Yukawa couplings the issue of 75 is likely to 
be more serious. For the Higgs boson self-coupling, on 
the other hand, one has to consider four-point functions 
and thus a mapping to massless two-point functions in- 
troduces infra-red divergences already at one-loop order. 
A promising method would be to assign a common mass 
to all fields and set all external momenta to zero (see, 
e.g., Ref. 0). 

Acknowledgements. We are grateful to Konstantin 
Chetyrkin and Johann Kiihn for useful comments, and 
we thank David Kunz for tests concerning the interface 
between FeynArts and q2e. This work is supported 
by DFG through SFB/TR 9 "Computational Particle 
Physics" . 



[1] 



D. J. Gross and F. Wilczek, Phys. Rev. Lett. 30 (1973) 
1343. 

H. D. Politzer, Phys. Rev. Lett. 30 (1973) 1346. 

D. R. T. Jones, Nucl. Phys. B 75, 531 (1974). 
W. E. Caswell, Phys. Rev. Lett. 33 (1974) 244. 

O. V. Tarasov and A. A. Vladimirov, Sov. J. Nucl. Phys. 

25 (1977) 585 [Yad. Fiz. 25 (1977) 1104]. 

E. Egorian and O. V. Tarasov, Teor. Mat. Fiz. 41 (1979) 

26 [Theor. Math. Phys. 41 (1979) 863]. 

O. V. Tarasov, A. A. Vladimirov and A. Y. Zharkov, 
Phys. Lett. B 93 (1980) 429. 

S. A. Larin and J. A. M. Ver maseren, Phys. Lett. B 303 
(1993) 334 jhep-ph/9302208] . 

T. van Ritbergen, J. A. M. Vermaseren and S. A. Larin, 
Phys. Lett. B 400 (1997) 379 [ hep-ph/970l"390] . 
Czakon, 



M 



485 



Nucl. Phys. B 710 (2005) 

Rev. D 59 (1999) 054005 



B 387 (1992) 373 [Erratum-ibid. 
|arXiv:hep-ph/0111190] . 



Phys. 
B 504 (1997) 551] 



[21 
[22 

[23 

[24 

[25 

[26 
[27 

[28 
[29 
[30 



arXi v:hep-ph/0411261] . 
M. Steinhauser, Phys 
[arXiv:hep-ph/9809507] . [31 
D. J. Gross and F. Wilczek, Phys. Rev. D 8 (1973) 3633. 
D. R. T. Jones, Phys. Rev. D 25 (1982) 581. [32 
M. S. Fischler and C. T. Hill, Nucl. Phys. B 193 (1981) 
53. 

M. E. Machacek and M. F. Vaughn, Nucl. Phys. B 222 [33 

(1983) 83. 

I. Jack and H. Osborn, Nucl. Phys. B 249 (1985) 472. [34 
M. E. Machacek and M. T. Vaughn, Nucl. Phys. B 236 

(1984) 221. 

M. Fischler and J. Oliensis, Phys. Lett. B 119 (1982) [35 
385. 

M. E. Machacek and M. F. Vaughn, Nucl. Phys. B 249 

(1985) 70. 

C. Ford, I. Jack and D. R. T. Jones, Nucl 



[36 



M. x. Luo and Y. Xiao, Phys. Rev. Lett. 90 (2003) 011601 



arXiv:hep-ph/0207271 
A. G. M. Pickering, J. A. Gracey and D. R. T. Jones, 
Phys. Lett. B 510 (2001) 347 [Ph ys. Lett. B 512 (2 001 
ERRAT,B535,377.2002) 230] [arXiv:hep-ph/0104247] . 
S. A. Larin, F. V. Tkachov and J. A. M. Vermaseren, 
preprint NIKHEF-H-91-18 (1991). 

M. E. Peskin and D. V. Schroeder, Reading, USA: 
Addison- Wesley (1995) 842 p 

L. N. Mihaila, J. Salomon, M. Steinhauser, in prepara- 
tion. 

P. Nogueira, J. Comput. Phys. 105 (1993) 279. 
R. Harlander, T. Seidensticker and M. Steinhauser, Phys. 
Lett. B 426 (1998) 125 [hep-ph/9712228] . 
T. Seidensticker, hep-ph/9905298 



J. A. M. Vermaseren, arXiv:math-ph/0010025 
T. Hahn, Comput. Phys. Commun. 140 (2001) 418 
[hep-ph/00i2260] . 

N. D. Christensen and C. Duhr, Comput. Phys. Com- 
mun. 180 (2009) 1614 |arXiv:0806.4l94l [hep-ph]]. 
V. A. Smirnov, "Applied asymptotic expansions in mo- 
menta and masses," Springer Tracts Mod. Phys. 177 
(2002) 1. 

M. Misiak and M. Miin z, Phys. Lett. B 344 (1995) 308 
[arXiv:hep-ph/9409454"] . 

A generalization of the final result to non-vanishing 
Yukawa couplings for the first and second generation is 
straightforward . 
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The two-point functions are infra-red safe. 



